REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  Information  Is  estimated  to  average  1  hour  per  response.  Including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  Information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  Information,  Including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson 
Davis  Highway,  Suite  1284,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503 


1.  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

9/19/97 


4.  TITLE  AND  SUBTITLE 

Waves  in  Plasma  Sheaths  and  at  Boundaries: 
Theory  and  Computer  Experiments 


6.  AUTHOR{S) 

Charles  Birdsall 


3.  REPORT  TYPE  AND  DATES  COVERED 

Final  -  7/1/94  -  9/30/97 


5.  FUNDING  NUMBERS 

N00014-94-1-1033 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS{ES) 
University  of  California,  Berkeley 
Electronics  Research  Laboratory 
253  Cory  Hall 
Berkeley,  CA  94720 


9.  SPNSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS<ES) 

Office  of  Naval  Resarch 

Charles  Roberson,  Scientific  Officer:  Code  312AI 
800  North  Quincy  Street 
Arlington,  VA  22217-5660 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

442427-23057 


10.  SPONSORING  /MONITORING 
AGENCY 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


OK  K 


j|£pzoTod  ie&  BuWic  mlooafil 


13.  ABSTRACT  (Maximum  200  words) 

The  goal  of  the  research  is  to  find  the  mechanisms  in  the  presheath,  sheath, 
and  wall  region  that  reflect,  transmit,  and  absorb  ion  acoustic  waves.  Ion  waves 
are  observed  to  propagate  from  the  bulk  plasma  into  and  completely  through  the 
sheath  to  the  wall. 


14.  SUBJECT  TERMS 

ion  acoustic  waves,  presheath,  sheath. 

wall  region. 

15.  NUMBER  OF  PAGES 

4 

crossed-f ield  diodes 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

18.  SECURITY  CLASSIRCATION 

19.  SECURITY  CLASSIFICATION 

20.  LIMITATION  OF  ABSTRACT 

OF  REPORT 

OF  THIS  PAGE 

OF  ABSTRACT 

unclassified 

unclassified 

unclassified 

NSN  7540-01-280-5500 


Standard  Form  298  (890104  Draft) 

Pr«*orlb«d  by  ANSI  Std.  239-18  298-01 


UNIVERSITY  OF  CALIFORNIA,  BERKELEY 


BERKELEY  •  DAVIS  •  IRVINE  •  LOS  ANGELES  •  RIVERSIDE  •  SAN  DIEGO  •  SAN  FRANCISCO 


SANTA  BARBARA  •  SANTA  CRLfZ 


COLLEGE  OF  ENGINEERING  BERKELEY,  CA  94720-1774 

ELECTRONIC  RESEARCH  LABORATORY 
253  CORY  HALL 
(510)  642-2301 

http://www.eecs.berkeley.edu/~erl 


Waves  in  Plasma  Sheaths  and  at  Boundaries: 
Theory  and  Computer  Experiments 


AASERT  Final  Report 

N00014-94- 1-1033 
7/1/94-9/30/97 


University  of  California,  Berkeley 
Electronics  Research  Laboratory 

Charles  K.  Birdsall,  Principal  Investigator 
Keith  Cartwright,  Graduate  Student  Researcher 


Office  of  Naval  Research 
9/30/97 
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Keith  Cartwright,  Graduate  Student  Researcher, 

Electronics  Research  Laboratory, 

University  of  California,  Berkeley 

This  is  a  description  of  the  work  Keith  Cartwright  has  done  as  a  graduate  student  in  the  last  year 
supported  by  an  AASERT  in  Professor  C.  K.  Birdsall’s  Plasma  Theory  and  Simulation  Group.  The 
four  main  areas  of  research  are:  (1)  the  reflection  of  ion  waves  from  the  sheath,  (2)  crossed-field  diodes, 
(3)  development  of  a  PIC  ion  -  Boltzmann  electrons  simulation  model,  and  (4)  the  development  of 
the  OOPIC  code. 

1.  The  goal  of  this  work  is  to  And  the  mechanisms  in  the  presheath,  sheath,  and  wall  region  that 
reflect,  transmit,  and  absorb  ion  acoustic  waves.  Ion  waves  are  observed  to  propagate  from  the 
bulk  plasma  into  and  completely  through  the  sheath  to  the  wall.  The  sheaths  used  for  study  have 
been  formed  next  to  a  metal  wall  and  around  a  grid.  These  waves  have  a  frequency  less  than 
the  ion  plasma  frequency.  Ion  waves  have  been  purposely  launched  as  well  as  spontaneously 
generated  by  an  ion  acoustic  instability.  A  continuous  spectrum  of  spontaneously  generated 
ion  acoustic  standing  waves  are  observed  and  eigenfunctions  have  been  measured  in  Id  PIC 
simulations. 

2.  In  cross  field  diode  simulations  we  are  contrasting  2d  results  with  Id  results  already  obtained 
in  planar^  and  cylindrical^  geometries.  This  investigation  indicates  the  transverse  dimension  is 
important  in  delaying  the  onset  of  virtual  cathode  oscillations  for  currents  above  the  limiting 
current.  The  work  so  far  has  been  limited  to  the  behavior  of  diodes  at  5  ^  B^ull  and  B  =  Q, 
The  delay  is  only  observed  in  the  magnetized  cases. 

3.  Kinetic  simulation  of  plasmas  in  which  equilibrium  occurs  over  ion  timescales  poses  a  computa¬ 
tional  challenge  due  to  the  disparate  timescales  of  the  electron  plasma  frequency  (~  10®),  the  ion 
plasma  frequency  10^),  ion  transit  frequency  10^),  and  the  ionization  frequency  10^). 
Hybrid  electrostatic  PIC  algorithms  are  presented  in  which  the  electrons  reach  thermodynamic 
equilibrium  with  the  ions  each  time  step.  The  collision  frequency,  iyrn(v),  can  be  a  tabulated  or 
fitted  function.  These  approximations  neglect  eflFects  faster  than  ion  time-scales,  decreasing  the 
computer  time  used  by  over  an  order  of  magnitude;  however,  they  increase  the  complexity  of  the 
boundary  conditions  and  the  simulation  is  no  longer  self-consistent.  Theoretical  ramifications 
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of  these  approximations  are  examined,  and  results  are  compared  with  full  PIC  simulations  of 
undriven  sheaths,  ion  acoustic  waves,  plasma  display  panels,  and  DC  discharges. 

4.  The  objectives  here  is  to  continue  development  the  OOPIC  PIC-MCC  simulation  program,  a 
particle-in-cell,  electromagnetic/electrostatic  code  which  models  devices  ranging  from  high  pres¬ 
sure  discharges  to  vacuum  electronics  to  microwave-beam  devices.  This  program  is  written  in 
Cd-+,  using  an  object-oriented  design  philosophy.  Most  of  this  work  has  been  done  in  collabo¬ 
ration  with  Dr,  John  Verboncoeur.  Our  accomplishments  are: 

(a)  Improvements  in  the  nonlinear  Boltzmann  electron  field  solver  in  OOPIC  code.  For  details 
see  Section  3.  It  is  important  to  extend  the  usefulness  of  OOPIC  from  high  frequency 
electromagnetic  models  and  electron  dominated  electrostatic  systems  to  the  domain  where 
low  frequency  ion  dynamics  are  dominant. 

(b)  A  Multigrid  Poisson  Solve  has  been  added. 

The  following  is  a  list  of  abstracts  of  presentations  made: 

1.  Refraction  and  Reflection  of  Ion  Acoustic  Solitons  by  Space  Charge  Sheath  formed  by  a 

Grid  * 

K.  L.  Cartwright,  and  C.  K.  Birdsall 

Experiments  have  shown^’"^  that  ion  acoustic  solitons  tunnel  through  the  space  charge  sheath  in  front  of  a  grid  without  time 
delay.  They  are  absorbed  resonantly  when  the  spatial  width  of  the  wave  is  close  to  the  characteristic  gradient  scale  length 
of  the  sheath.  The  reflection  and  transmission  coefficients  found  in  these  experiments  have  compared  well  with  theory  in 
the  long  wavelength  limit.®  However,  to  achieve  this  comparison,  two  parameters  were  added  that  were  not  in  the  original 
theory.  These  parameters  allow  for  the  absorption  of  wave  energy  by  the  space  charge  sheath.  The  goal  of  our  numerical 
simulations,  designed  to  reproduce  the  experimental  results,  is  to  uncover  the  mechanism  of  this  energy  loss  and  large  speed 
of  propagation  through  the  sheath. 

2.  Hybrid  PIC  Acceleration  Schemes  for  Discharges  ^ 

K.  L.  Cartwright,  J.  P.  Verboncoeur,  and  C.  K.  Birdsall 

Kinetic  simulation  of  plasmas  in  which  equilibrium  occurs  over  ion  timescales  poses  a  computational  challenge  due  to  the 
disparate  timescales  of  the  electron  plasma  frequency  10^),  the  ion  plasma  frequency  10^),  ion  transit  frequency 
10®),  and  the  ionization  frequency  10^).  Hybrid  electrostatic  PIC  algorithms  are  presented  in  which  the  electrons 
reach  thermodynamic  equilibrium  with  the  ions  each  time  step.  There  are  three  different  approximations  for  the  electrons. 
First,  the  nonlinear  Boltzmann  relationship  for  the  electrons  can  be  applied  to  the  bulk  of  a  plasma.  Second,  there  is  a 
truncated  Maxwellian  which  is  used  in  undriven  sheaths;  this  approximation  truncates  the  electron  distribution  at  the  wall 
potential.  Progress  is  shown  for  a  method  that  will  be  used  in  a  DC  electric  field.  This  method  aissumes  a  steady-state 
electron  distribution  function  in  a  DC  electric  field,  with  elastic  collisions  between  electrons  and  neutral  gas  atoms.  The 
collision  frequency,  can  be  a  tabulated  or  fitted  function;  the  method  is  implemented  with  a  constant  cross  section 

*  American  Physical  Society  Division  of  Plasma  Physics,  November  1995,  Denver,  Co. 

^lEEE  International  Conference  on  Plasma  Science,  May  1997,  San  Diego,  CA. 
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and  is  planned  to  be  implemented  using  Argon  cross-sections.  These  approximations  neglect  effects  faster  than  ion  time- 
scales,  decreasing  the  the  computer  time  used  by  over  an  order  of  magnitude;  however,  they  increase  the  complexity  of 
the  boundary  conditions  and  the  simulation  is  no  longer  self-consistent.  Theoretical  ramifications  of  these  approximations 
are  examined,  and  results  are  compared  with  full  PIC  simulations  of  undriven  sheaths  and  preliminary  comparisons  of  ion 
acoustic  waves  and  DC  discharges. 

3.  Refraction  and  Reflection  of  Ion  Acoustic  Solitons  by  Space  Charge  Sheath  ^ 

K.  L.  Cartwright,  and  C.  K.  Birdsall 

Experiments  have  shown, that  ion  acoustic  solitons  tunnel  through  the  space  charge  sheath  in  front  of  a  grid  without  time 
delay.  They  are  absorbed  resonantly  when  the  spatial  width  of  the  wave  is  close  to  the  characteristic  gradient  scale  length 
of  the  sheath.  The  reflection  and  transmission  coefficients  found  in  these  experiments  have  compared  well  with  theory  in 
the  long  wavelength  limit.®  However,  to  achieve  this  comparison,  two  parameters  were  added  that  were  not  in  the  original 
theory.  These  parameters  allow  for  the  absorption  of  wave  energy  by  the  space  charge  sheath.  The  goal  of  our  numerical 
simulations,  designed  to  reproduce  the  experimental  results,  is  to  uncover  the  mechanism  of  this  energy  loss  and  large  speed 
of  propagation  through  the  sheath. 

4.  Relaxation  of  Virtual  Cathode  Oscillations  due  to  Transverse  Effects  in  a  Crossed-Field 

Diode  ^ 

K.  L.  Cartwright,  J.  P.  Verboncoeur,  V.  P.  Gopinath,  and  C.  K,  Birdsall 

Recent  studies  of  cylindrical  and  planar^  crossed-field  diodes  indicate  the  transverse  dimension  plays  a  role  in  delaying 
the  onset  of  virtual  cathode  oscillations  for  currents  injected  above  the  theoretical  one-dimensional  limiting  current.  For  Id 
and  2d  planar  devices,  the  limiting  current  for  the  magnetized  (R  =  and  unmagnetized  diodes  is  examined  for  cold 

and  thermal  injection.  A  significant  difference  between  the  Id  and  2d  smooth  wall  diodes  is  that  the  transverse  direction 
provides  an  extra  degree  of  freedom  which  allows  the  electrons  to  warm  rapidly.  The  mechanism  of  this  warming  appears 
to  be  an  instability  in  the  transverse  direction.  The  simulations  show  three  different  ‘states’;  laminar  flow,  virtual  cathode 
oscillation  and  warm  flow.  Warm  flow  occurs  when  the  electrons  have  a  spread  of  energy,  either  due  to  an  instability  or 
by  thermal  injection,  when  they  pass  though  the  potential  minimum.  Birdsall  and  Bridges®  showed  that  a  small  thermal 
spread  (few  percent)  of  injected  electrons  damps  virtual  cathode  oscillations  in  a  short  circuited  diode.  This  warming  effect 
allows  ‘warm  flow’  to  exist  on  the  2d  state  diagram  {J/Jc  versus  B/Bhuu)  which  is  not  found  on  the  Id  state  diagram  for 
cold  emission.  Parameter  space  is  explored  on  these  state  diagrams  for  B  =  0  and  B  =  BhuU  for  J  near  state  transitions 
{.J  ~  these  PIC  simulations,  the  particles  are  injected  in  such  a  way  as  to  keep  the  transverse  density  as  uniform 

as  possible.  A  quiet  particle  injection  algorithm  was  used  so  that  particle  noise  does  not  warm  in  the  transverse  direction. 
Time  centering  of  the  velocities  of  the  injected  particles  was  also  examined. 


^Sept  1996,  Sendai  Japan. 
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